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Abstract. Since the first limit on the (local) primordial non-Gaussianity parameter, 

. /nl, was obtained from the Cosmic Background Explorer (COBE) data in 2002, 

| observations of the cosmic microwave background (CMB) have been playing a central 

' role in constraining the amplitudes of various forms of non-Gaussianity in primordial 

O !■ fluctuations. The current 68% limit from the 7-year data of the Wilkinson Microwave 



Anisotropy Probe (WMAP) is /nl = 32 ± 21, and the Planck satellite is expected to 



!— i ■ reduce the uncertainty by a factor of four in a few years from now. If /nl 3> 1 is found 

by Planck with high statistical significance, all single-field models of inflation would 
be ruled out. Moreover, if the Planck satellite finds /nl ~ 30, then it would be able 
to test a broad class of multi-field models using the four-point function (trispectrum) 
\ test of tnl > (6/nl/5) 2 . In this article, we review the methods (optimal estimator), 

j^. ■ results (WMAP 7-year), and challenges (secondary anisotropy, second-order effect, and 

, foreground) of measuring primordial non-Gaussianity from the CMB data, present a 

science case for the trispectrum, and conclude with future prospects. 
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1. Introduction 

The physics of the very early, primordial universe is best probed by measurements of 
statistical properties of primordial fluctuations. The primordial fluctuations are the 
seeds for the temperature and polarization anisotropics of the CMB and the large-scale 
structure of the universe that we observe today. Therefore, both the CMB and the 
large-scale structure are excellent probes of the primordial fluctuations. In this article, 
we shall focus on the CMB. See the article by V. Desjacques and U. Seljak in this volume 
for the corresponding review on the large-scale structure as a probe of the primordial 
fluctuations. 

This article reviews a recent progress on our using the CMB as a probe of 
a particular statistical aspect of primordial fluctuations called u non-Gaussianity." 
Reviews on this subject were written in 2001 [1] and 2004 J3|. The former review 
would be most useful for those who are new to this subject. 

In this article, we focus on the new discoveries that have been made since 2004. 
Particularly notable ones include: 
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1. It has been proven that a//inflation models (not just simple ones jl, Q) based upon 
a single scalar field would be ruled out regardless of the details of models j^], if the 
primordial non-Gaussianity parameter called /nl (more precisely, the "local type" 
/nl as described later) is found to be much greater than unity. 

2. The optimal method for extracting /nl from the CMB data has been developed 
0-0] and implemented 10]. The latest limit on the local-type /nl from the WMAP 
7- year temperature data is /nl = 32 ± 21 (68% CL) [ill ]. 

3. The most serious contamination of the local- type /nl due to the secondary CMB 
anisotropy, the coupling between the Integrated Sachs- Wolfe (ISW) effect and the 
weak gravitational lensing, has been identified jil, ll2T-[l6l] . However, note that the 
astrophysical contamination such as the Galactic foreground emission and radio 
point sources may still be the most serious contaminant of /nl- These effects 
would pose a serious analysis challenge to measuring / N l from the Planck data. 

4. The importance of distinguishing different triangle configurations of the three-point 
function of the CMB was realized 17J, ll8j| and has been fully appreciated. It 
has been shown by many researchers that different configurations probe distinctly 
different aspects of the physics of the primordial universe. The list of possibilities 
is long, and a terribly incomplete list of references on recent work (since ~ 
2004) is: [19| on a general analysis of various shapes; [20l-|25| on the local shape 



(k 3 <C ki m k 2 ); |l7l. I26M28I] on the equilateral shape (ki m k 2 ~ k 3 ); [29|, |30j on the 



flattend (or folded) shape {k\ w 2k 2 ~ 2/c 3 ); [3lJ,|32| on the orthogonal shape (which 



33 



-37 on combinations 



is nearly orthogonal to both local and equilateral shapes); 
of different shapes; and 38|, |39| on oscillating bispectra. Also see references therein. 

The connected four-point function of primordial fluctuations has been shown to be 
an equally powerful probe of the physics of the primordial universe. In particular, 
a combination of the three- and four-point functions may allow us to further 
distinguish different scenarios. Many papers have been written on this subject 



over the last few years: |40|-|43j on single-field models; |35|, |44j-|69[ on multi-field 
models; and [70, |7l| on isocurvature perturbations. CMB data are expected to 
provide useful limits on the parameters of the "local-form trispectrum," r NL and 
<?nl 72|, |73| . Preliminary limits on these parameters have been obtained from the 
WMAP data by 0,0- 

The number of researchers working on primordial non-Gaussianity has increased 
dramatically: Science White Paper on non-Gaussianity submitted to Decadal Survey 
Astro2010 was co-signed by 61 scientists 76 . 



2. Gaussian versus non-Gaussian CMB anisotropy 

2.1. What do we mean by "Gaussianity"? 

What do we mean by "Gaussian fluctuations"? Let us consider the distribution 
of temperature anisotropy of the CMB that we observe on the sky, AT(n). The 
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temperature anisotropy is Gaussian when its probability density function (PDF) is given 
by 



P(AT) 



1 



(2 7 r) Ar p^/ 2 |e| 1/2 



exp 



(1) 



where ATj = AT(n), ^ = (ATjATj) is the covariance matrix (or the two-point 
correlation function) of the temperature anisotropy, |£| is the determinant of the 
covariance matrix, and iVpi x is the number of pixels on the sky. 

We often work in harmonic space by expanding AT using spherical harmonics: 
AT(n) = Y J im a imYim{n). The PDF for a Lm is given by 



P(a) 



1 



(2 7r )AWm/2|(7|l/2 



exp 
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where Ci m y m i = (a>i m (ii>m') , and iV harm is the number of / and m. When a\ m is statistically 
homogeneous and isotropic (which is not always the case because of, e.g., non- uniform 



noise), one finds C, 



Im.l'm' 



P(a) 



n 

i in 



Cid~u>5 mrn i } and thus the PDF simplifies to 

-\a lm \y(2d) 



(3) 



Here, Ci is the angular power spectrum. The latest determination of C\ of the CMB 
temperature anisotropy is shown in Figure [TJ 

The important property of a Gaussian distribution is that the PDF is fully specified 
by the covariance matrix. In other words, the covariance matrix contains all the 
information on statistical properties of Gaussian fluctuations. When the PDF is given 
by equation [3j the power spectrum, C[, contains all the information on a/ m . This is 
not true for non-Gaussian fluctuations, for which one needs information on higher-order 
correlation functions. 

Let us close this subsection by noting that a non-zero deviation of the covariance 
matrix from the diagonal form, AC/ mj // m / = C/ mj // m / — Ci5w5 mm ', does not imply non- 
Gaussianity: the PDF can be a Gaussian with a non-diagonal covariance matrix as 
given in equation |2j A non-zero ACi m ^ m i may arise in cosmological models that violate 
statistical isotropy. Such models may yield anisotropic Gaussian fluctuations; thus, one 
must distinguish between non-Gaussianity and a violation of statistical isotropy. 



2.2. What do we mean by "non-Gaussianity" ? 

What do we mean by "non-Gaussian fluctuations"? Any deviation from a Gaussian 
distribution (such as equation [[] or |2]) is called non-Gaussianity. When fluctuations 
in the CMB are non-Gaussian, one cannot generally write down its PDF, unless one 
considers certain models (e.g., inflation). Nevertheless, when non-Gaussianity is weak, 



one may expand the PDF around a Gaussian distribution 



80 
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Figure 1. The angular power spectrum of the CMB temperature anisotropy, Cz, 
measured from the WMAP 7- year data 77| , along with the temperature power spectra 
from the ACBAR [3] and QUaD [l^] experiments. The solid line shows the best-fitting 
6-parameter flat ACDM model to the WMAP data alone. The angular power spectrum 
contains all the information on fluctuations in the CMB, if fluctuations are Gaussian. 
If fluctuations are non-Gaussian, one must use the higher-order correlation functions 
(such as three- and four-point functions) to fully ex ploi t the cosmological information 
contained in the CMB. This figure is adopted from [111 ]. 
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2E|mEl'm' a *m( C 1 )l m ,l'm' a l'm' 



(2 7r )iV hami /2| ( 7|l/2 



(4) 



Here, the expansion is truncated at the three-point function (bispectrum) of a lm , 
and thus we have assumed that the connected four-point and higher-order correlation 
functions are negligible compared to the power spectrum and bispectrum. (This 
condition is not always satisfied.) By evaluating the above derivatives, one obtain^""" 



P(a) 



(27r) JV '»>«™/ 2 |C7| 1 /i 



exp 



_ 2 ^ ^ °* lm ^ 1 )lm,l'm'0,Vm' 
Ira I'm 1 



x < 1 + g ( a hmi a l2m2 a l3m 3 ) (C a)l imi (C 1 a)/ 2m2 (C 1 o)i 3 



<H3 



all liirij 

— 3(C 1 )i 1 m 1 ,i 2 m 2 {C 1 a)i 3m , i |. (5) 

This formula is useful, as it tells us how to estimate the angular bispectrum, 
( a i 1 mi a >hm2 a i3m 3 ) , optimally from a given data by maximizing this PDF. In practice, 
we usually parametrize the bispectrum using a few parameters (e.g., /nl), and estimate 
those parameters from the data by maximizing the PDF with respect to the parameters. 

| Babich (8lj | derived this formula for Ci m i' m i = CidwS mrn i . 
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In the limit that the contribution of the connected four-point function (trispectrum) 
to the PDF is negligible compared to those of the power spectrum and bispectrum, 
equation [5] contains all the information on non-Gaussian fluctuations characterized 
by the covariance matrix, Ci imi> i 2m2 = (az* imi az 2 m 2 ) > an d the angular bispectrum, 
m l a 'i2m2 a i3,ms) ■ This approach can be extended straightforwardly to the trispectrum 
if necessary. 



3. Extracting /nl from the CMB data 

3.1. General formula 

We have not defined what we mean by "/nl-" For the moment, let us loosely define it 
as the amplitude of a certain shape of the angular bispectrum: 

i 

where the function, &^] 2 ; 3 , is called the "reduced angular bispectrum" [82] and defines 
the shape of the angular bispectrum for a given model denoted by an index i (which 
may refer to, e.g., "local," "equilateral," "orthogonal," etc), and G™^™ 3 ls ^ ne so-called 
Gaunt integral, defined by 

The physical role of the Gaunt integral is to assure that (h,mi), (h,^), and (7 3 ,m 3 ) 
form a triangle. In the small- angle limit, the Gaunt integral becomes a 2-d delta 



function: G^Z^ -> (2vr) 2 ^(Z 1 + Z 2 + Z 3 ) [83| 



Given this parametrization, one can then maximize the PDF given in equation [5] 
with respect to by solving d In Pj df^l = and find the optimal estimator: 

& = EC* 1 " V*- (s) 

j 

One can use this formula to determine multiple amplitudes of angular bispectra 
simultaneously. 

Here, Si are given by the data as 



1 

all Im 



C — x \ ~" r'rn 1 m 2 m 3 i(i) 

y hhh °hhh 



x f(C 1 a)i imi (C 1 a)i 2Tn2 (C 1 a)i zm . i — 3(C 1 )i imi ,i 2m2 (C 1 a)i sm . i , (9) 

where 1/6 is included such that in equation [8] becomes the Fisher matrix of fjQ. In 
other words, the covariance matrix of /^l is given by the inverse of F^, i.e., 

(F-% = (fSif}ii)- (fSb (Ai). (io) 

The 68% uncertainty in /^l i s gi ven by A/^ = (F _1 )jj. 

Using the definition of the Gaunt integral given in equation [3 we rewrite Si as 
1 ( 

Si = 6 J d2fl ^ 6 ^ 2 ' 3 [ eZl ^) e ^(^) e ^(^) - 3d W2 (n)e, 3 (n)] , (11) 

hhh 
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where 

e,(n) = ^(C^a^mYUh), (12) 

m 

dw(n) = Y.^'^i^i'm'Yirni^Yi^^n). (13) 

mm' 

Here, the summation over m can be done using the Fast Fourier Transform (FFT), as 
Yi m (0,(f>) oc e tm ^. This technique is used by the HEALPix package [841 ]. and thus one 
may use HEALPix to do this summation. To compute du'(n), one may use Monte Carlo 
simulations. Namely, as C/ m ^/ m / = (a* m aj' m '), we have the exact relation between du> 
and e/: dw(h) = (ei(h)ei>(h)). One can evaluate the ensemble average using the Monte 
Carlo simulation of the CMB and the instrumental noise. Let us denote this operation 
by div (n) = {Gi(n)eii(h)) MC . The final formula for Si is 
1 ( 

Si = 6 J d2fl ^ h ^ hh M™K(n)e/ 3 (n) - 3e h {h)(e h (h)e h (h)) MC } (14) 

which is valid for general forms of &[*l 2 j 3 - Note that the integral, / d 2 n, must be done over 
the full sky, even in the presence of the mask: the information on the mask is included 
in the calculation of the Fisher matrix, F^. The only assumptions that we have made 
so far are: (1) each angular bispectrum component has only one free parameter, i.e., 
the amplitude, and (2) non-Gaussianity (if any) is weak, and the PDF of ai m is given 
by equation [51 

Finally, the explicit form of the Fisher matrix is given by 

771 _ /sky \ ^ \ -» r>m 1 m 2 m 3 Ai) 
all Im all I'm' 

x (C 1 )l 1 m 1 ,l' 1 m' 1 {C 1 )l 2 m 2 ,l' 2 m' 2 {C *) 7 3 m 3 'l^l'JL'^ 3 ' (^) 

where / s k y is the fraction of the sky outside of the mask. When the covariance matrix 
is diagonal, the expression simplifies to 

',(*) uU) 



r LW 1X1 J 

rp _ isky \ - T u hl 2 h u hhh 

b all I ^h^h^te 



where 



T = (rrnim2mz\2 _ (2/j + l)(2/ 2 + l)(2/ 3 + 1) / h h h 

lhhh ~ [ ^ ' ~ 4tt I 

all m \ 




Equation [16] may also be written as 



6 « 



^ij - /sky 2^ L hhh n n n \ ' 

h<h<h ^h^h^h^hhh 

where A^^jg = 1, 2, and 6 when all of /j's are different, two of l[ are the same, and all 
of l^s are the same, respectively. 
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3.2. Poisson bispectrum 



As a warm up, let us consider the simplest example: point sources randomly distributed 
over the sky. As mentioned already, this is a contamination of the primordial non- 
Gaussianity parameters, and thus an accurate measurement of this component is 
quite important, especially for Planck as well as for the South Pole Telescope (SPT) 
and Atacama Cosmology Telescope (ACT), which are working at high frequencies 
[y > 100 GHz) where star-forming galaxies dominate a/ TO at I > 1000. 

For the Poisson distribution, the reduced bispectrum is independent of multipoles, 
bfj la = 1, in the absence of window functions, and is given by 

(19) 



rsrc 
°hhh 



w h w h wi 3 , 

in the presence of window functions. Here, W\ is an experimental window function (a 
product of the beam transfer function and the pixel window function). Let us then 
use 6 src (instead of /nl because this component has nothing to do with primordial 
fluctuations) to denote the amplitude of the Poisson bispectrum. 
From the data, we measure 5 src given by 
1 



S S r C = lj d2fl F 3 (^) " 3£(n)(£ 2 (n)> MC J , (20) 
where a map E(fi) is defined by j85| 

E {™) = X^e z (n) = ^ w;(C" 1 a) im Yi m (n), (21) 

( lm 

where e;(n) is given by equation (fT2l . The E(h) map is a Wiener-filtered map of point 
sources randomly distributed on the sky. 



3.3. Primordial bispectra 



(Most of this subsection is adopted from Section 6.1 of [111].) During the period of 
cosmic inflation |87H92| , quantum fluctuations were gen erated an d bec ame the seeds for 



93-961. See 97-1021 for reviews 



the cosmic structures that we observe today 

Inflation predicts that the statistical distribution of primordial fluctuations is nearly 
a Gaussian distribution with random phases. Measuring deviations from a Gaussian 
distribution, i.e., non-Gaussian correlations in primordial fluctuations, is a powerful 
test of inflation, as how precisely the distribution is (non-) Gaussian depends on the 



detailed physics of inflation. See [2|, |76] for reviews 



The observed angular bispectrum is related to the 3-dimensional bispectrum of 
primordial curvature perturbations, (CfciCfc 2 Cfc 3 ) = (Z^S (ki + k 2 + k 3 )B^(k 1 , k 2 , k 3 ). In 
the linear orde r, th e primordial curvature perturbation is related to Bardeen's curvature 

The CMB 

-k 



5$ 
3 V 



perturbation [1031 ] in the matter-dominated er a, by £ 
temperature anisotropy in the Sachs- Wolfe limit 1051 ] is given by AT/T 
We write the bispectrum of $ as 

($(fc 1 )$(fc 2 )$(fc 3 )) = (2 7 r) 3 5 D (fc 1 + k 2 + k 3 )F(h, k 2 , k 3 ). 



104]. 



.1$ 

3^ 



(22) 
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(a) squeezed triangle 
(k ~k >>k ) 

v 1 2 3' 



(b) elongated triangle 
(k =k +k ) 

V 1 2 3' 



(c) folded triangle 
(k =2k =2k ) 

V 1 2 3' 






(d) isosceles triangle 
(k >k =k ) 

V 1 2 3' 




(e) equilateral triangle 
(k =k =k ) 

v 1 2 3' 





Figure 2. Visual representations of triangles forming the bispectrum, £>$(fci, &2, £3), 
with various combinations of wavenumbers satisfying &3 < k% < A»i. This figure is 
adopted from [86| . 



There is a useful way of visualizing the shape dependence of the bispectrum. 
We can study the structure of the bispectrum by plotting the magnitude of 
F(ki, k 2 , k 3 )(k 2 /ki) 2 (k 3 /ki) 2 as a function of k 2 jk\ and ks/ki for a given k\, with 
a condition that ki > k 2 > k 3 is satisfied. In order to classify various shapes of 
the triangles, let us use the following names: squeezed {k\ ~ k 2 ^> k 3 ), elongated 
(k\ = k 2 + ^3), folded [k\ = 2k 2 = 2/c 3 ), isosceles {k 2 = fc 3 ), and equilateral 
(ki = k 2 = k 3 ). See (a)-(e) of Figure. [2] for the visual representations of these triangles. 

We shall explore 3 different shapes of the primordial bispectrum: "local," 
"equilateral," and "orthogonal." They are defined as follows: 

1. Local form. The local form bispectrum is given by 

-flocal(^l5 k 2 ,k 3 ) 
flocalfo I h. \p 



82 



106 



107 



2A 2 /n°l c 



1 



7 4-n s 7,4- 
ft 2 



+ P*(k 2 )P*(k 3 ) + P*(fc 3 )P*(fcl)] 

- + (2 perm.} 



(23) 



where P$ = A/fc 4- ™" is the power spectrum of $ with a normalization factor A. 
This form is called the local form, as this bispectrum can arise from the curvature 
perturbation in the form of $ = + /nl ^*i> where both sides are evaluated 



at the same location in space ($£ is a linear Gaussian fluctuation) |§| The local 
form, Fi oca i(ki,k 2 ,k 3 )(k 2 /ki) 2 (k 3 /ki) 2 , peaks at the so-called "squeezed" triangle 

§ However, $ = <J>l + /nl^^l is n °t the oru y wa Y to produce this type of bispectrum. One can also 
produce this form from multi-scalar field inflation models where scalar field fluctuations are nearly 
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Figure 3. Shapes of the primordial bispectra. Each panel shows the normalized 
amplitude of F(ki, k 2 , k 3 )(k2/k 1 ) 2 (k 3 /ki) 2 as a function oik 2 /ki and fc 3 /fci for a given 
k\, with a condition that k 3 < k 2 < k\ is satisfied. As the primordial bispectra 
shown here are (nearly) scale invariant, the shapes look similar regardless of the 
values of k\. The amplitude is normalized such that it is unity at the point where 
F(k\, k 2 , k 3 )(k 2 /ki) 2 {k 3 /ki) 2 takes on the maximum value. (Top Left) The local form 
given in equation 1231 which peaks at the squeezed configuration. Note that the most 
squeezed configuration shown here has k\ = k 2 = 100^3 . (Top Right) The orthogonal 
form given in equation 1271 which has a positive peak at the equilateral configuration, 
and a negative valley along the elongated configurations. (Bottom Left) The equilateral 
form given in equation 1261 which peaks at the equilateral configuration. Note that all 
of these shapes are nearly orthogonal to each other. 



for which k 3 <C k 2 
obtain 



ki [18|. See the top- left panel of Figure El In this limit, we 



(24) 



How large is /^l 8,1 from inflation? The earlier calculations showed that /^l 8,1 from 



single-field slo w-roll inflation would be of order the slow-roll parameter, e ~ 10 



106 



114 



1 15l | . More recently, Maldacena (^J and Acquaviva et al. [4] found that 
the coefficient of P$(/ci)P$(/c3) from the simplest single-field slow-roll inflation with 
the canonical kinetic term in the squeezed limit is given by 



-FlocalOl, h, k 3 -> 0) = -(1 



n> 



(25) 



Comparing this result with the form predicted by the /^n cal model, one obtains 



local 
NL 



(5/12)(l - n s ), which gives /^ al = 0.015 for n 



'NL 

0.963. 



scale invariant [22|; multi-scalar models called "curvaton" s cenarios [2dl . LLQ8(; multi- field models in 
which one field modulates the decay rate of infiaton field 21, 109l llo| : multi- field models in which a 



violent product ion of pa rticles and non-linear reheating, called "preheating," occur due to parametric 
resonances models in which the universe contracts first and then bounces 231. 
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2. Equilateral form. The equilateral form bispectrum is given by 
-Pequii(^i) k 2 ,k 3 ) 
= 6A 2 fZ ai 1 1 1 



-n s h.i-n 3 
l%2 



hv 2 



^3 



(M 2 ^3) 2{4 - ns)/3 



+ 



1 



,(4-n s )/3,2(4-n s )/3,4- 
rv 2 - 



h (5 perm. 



(26) 



.fi, 1 i% 2 ft-3 

This function approximates the bispectrum forms that arise from a class of inflation 
models in which scalar fields have non-canonical kine tic terms. One example is the 
so-called Dirac-Born-Infeld (DBI) inflation 26_|, ll 16 ] . which gives /nl" 1 ^ ~~ V c s 
in the limit of c s <C 1, where c s is the effective sound speed at which scalar 
field fluctuations propagate relative to^ the spe ed of lig ht. There are various 

The equilateral form, 



other models that can produce /nl" 1 19_|, |27J, 



28 



117 



118 



F eq uii(ki, k 2 , k 3 )(k 2 /ki) 2 (k 3 /ki) 2 , peaks at the equilateral configuration for which 
k\ = k 2 = k 3 . See the bottom- left panel of Figure El The local and equilateral 
forms are nearly orthogonal to each other, which means that both can be measured 
nearly independently. 

3. Orthogonal form. The orthogonal form, which is constructed such that it is 



nearly orthogonal to both the local and equilateral forms, is given by 31 



-^orthog 

(ki,k 2 ,k 3 ) 



2 i-orthog 



NL 



4— n» !_4— n. 



At 2 ft-3 



ft 3 At 1 



+ 



(4-n s )/3 fc 2(4-n s )/3,^4 



U4-n a 
fi, 3 



(5 perm. 



(27) 



This form approximates the forms that arise from a linear combination of 
higher-derivative scalar-field interaction terms, each of which yields forms 
similar to the equilateral shape. Senatore, Smith and Zaldarriaga [3l| found 
that, using the "effective field theory of inflation" approach 117J, a certain 
linear combination of similarly equilateral shapes can yield a distinct shape 
which is orthogonal to both the local and equilateral forms. The orthogonal 
form, F orthog (fc 1 , k 2 , k 3 ) (k 2 / ki) 2 (k 3 / ki) 2 , has a positive peak at the equilateral 
configuration, and a negative valley along the elongated configurations. See the 
top- right panel of Figure [31 

Note that these are not the most general forms one can write down, and there 
are other forms whic h would probe different aspects of the physics of inflation 



19 



29 



33 



37 



69 



119 



Of these forms, the local form bispectrum has special significance. Creminelli and 
Zaldarriaga [5] showed that not only models with the canonical kinetic term, but all 
single- inflation models predict the bispectrum in the squeezed limit given by equation 1251 
regardless of the form of potential, kinetic term, slow-roll, or initial vacuum state. Also 
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sec 



19l |28| . Ill7j . This means that a convincing detection of /nl would rule out all 



single-field inflation models. 



3.4- Optimal estimator for f^ al 

Given the form of $, one can calculate the harmonic coefficients of temperature and 
E'-mode polarization anisotropics as 

d 3 k 



a lr, 



Im 



47r(-z)' 



4tt(-z)' 



(27T) 



Mk)g Tl (k)Yr m (k), 



(l + 2)\ f d 3 k 



Mk)9Pi(k)Y lm (k), 



(2c 



(29) 



\(l-2)\J (2vr)3 

where gn{k) and gpi(k) are the radiation transfer functions of the temperature and 
polarization anisotropies, respectively, which can be calculated by solving the linearized 
Boltzmann equ ations. On e m ay use the publicly-available Boltzmann codes such as 
CMBFAST 120] or CAMB |l21 | for computing the radiation transfer functionsjjj] 

From now on, we shall focus on the temperature anisotropy, largely for simplicity. 
(See 122l . Il23| for the treatment of polarization in the angular bispectrum.) The limits 



on /nl expected from Planck are dominated by the temperature information, and thus 
the polarization information is not expected to yield competitive limits over the next, 
say, > 5 years. 

For the local-form bispectrum given in equation [231 the reduced bispectrum 
(equation [6]) is given by [82 

= 2 / r " dr iPh(r)P h (r)a h (r) + (2 perm.)] w h w h w h , (30) 

where 

2 



2 f 

a i( r ) = ~ k 2 dk 9Ti(k)ji(kr) 

7T J 



A(r) 



7T 



k 2 dk P*(k)gn{k)ji{kr). 



(31) 
(32) 



Using this form in equation [TU one finds <Si oca i as 



5, 



local 



r dr / a I n 



A(h,r)B 2 (h,r) - 2B(h)(A(n,r)B(h,r)) 



MC 



-A(h, r)(B 2 (h, r)) 



MC 



(33) 



which can be measured from the data. Here, maps A(n,r) and B(n,r) are defined by 



A(h,r) = ^wiai(r)ei(h) = ^wiai(r)(C 1 a)i m Yi m (h) 



l in 



B(h,r) = ^wif3i(r)ei(h) = J2 w iPi( r )( c 1 a)i m Y lm (n), 



(34) 
(35) 



I'm 



where e^(n) is given by equation ffl2|) . 

|| A CMBFAST-based code for computing gri(k) and gpi(k) is available at 
http://gyudon.as.utexas.edu/~komatsu/CRL, A recent version of CAMB has an option to cal- 
culate these functions ( |http: / /camb.info| ). 
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3. 5. Optimal estimator for f^ d 

For the equilateral-form bispectrum given in equation [26], the reduced bispectrum is 
given by [71] 

€S 3 = - 3 ftl + 6 / r * dr [MrHAr)S h (r) + (5 perm.) 

-25 h (r)5 h (r)5 h (r)] w h w h w hl (36) 

where 



7t(r) = - / k 2 dk Pl /3 (k)g Tl (k) 3l (kr) } (37) 

7T J 

Sl (r) = - / fc 2 d/c P^i^gnik^ikr). (38) 
Using this form in equation [TJJ one finds Sequii as 

^e q uii= -3S loca i + 6 J r 2 dr J d 2 h {B(n,r)C(n,r)D(n,r) 

- B{fi)(C{n,r)D{h,r)) M c - C{n)(B{n,r)D{h,r)) MC 
-D(n){B{h, r)C(n, r)) MC - ^ [^ 3 (n, r) - 3£>(n, r)(Z> 2 (n> r)) MC ] } 



(39) 



which can be measured from the data. Here, maps C(n,r) and D(h,r) are defined by 


C(n,r) = ^u7 / 7;(r)e i (n) = ^ wai(r){C' l a) lm Yi m (n), (40) 

l Im 

D(n } r) = ^ Wi Mr)e ( (n) = ^(^(C^aV^n). (41) 



5. 5. Optimal estimator for /nl ° g 

For the equilateral-form bispectrum given in equation [271 the reduced bispectrum is 
given by [31] 



CP = 36Effi - 12 J r " dr ^r)K{r)Wv h w h . (42) 
Using this form in equation [TU one finds .Sorthog as 

Orthog = 3S cquil - 2 J r 2 dr J d 2 h [D 3 (h,r) - 3D(h, r)(D 2 (h, r)) MC ] ,(43) 
which can be measured from the data. 

4. Secondary anisotropy 

^.1. General formula for the lensing- secondary coupling 

Given the special importance of the local-form bispectrum, we must understand what 
other (non-primordial) effects might also produce the local form, potentially preventing 
us from measuring /nl ■ 
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The local-form bispectrum is generated when the power spectrum of short- 
wavelength fluctuations is modulated by long-wavelength fluctuations; thus, a 
mechanism that couples small scales to large scales can potentially generate the local- 
form bispectrum. 

The weak gravitational lensing provides one such mechanism. The local-form 
bispectrum may then be generated when long- and short-wavelength fluctuations are 
coupled by the lensing. To see how this might happen, let us write the observed 
temperature anisotropy in terms of the original (unlensed) contribution from the last 
scattering surface at z = 1090, AT P (where "P" stands for "primary" 1 ), the lensing 
potential, 0, and the secondary anisotropy generated between z = 1090 and z — 0, AT S 
(where "5"' stands for " secondary"): 

AT(n) = AT p (h + d<j>) + AT s (h) 

« AT p (h) + \{d<P) ■ (dAT p )](h) + AT s (h), (44) 



where 



in 



-2 



dr- 



rr % 



-<3>(r, hr) 



(45) 



and r* is the comoving distance out to z — 1090, and $ is Bardeen's curvature 
perturbation, which is related to the usual Newtonian gravitational potential by 
$ = — $Ncwton- Transforming this into harmonic space and computing the reduced 
bispectrum, one obtains [121 ] 



Liens— S 
°hhh 



Mh + i)-faft + i) + W. + i) g P Cf + (5 perm , 



w h w h w h , (46) 



where Cf is the power spectrum of the CMB from the decoupling epoch only (i.e., no 
ISW), and Cf s = (4>J m af m ) is the lensing-secondary cross-correlation power spectrum. 

From this result, one finds that a non-zero bispectrum is generated when the 
secondary anisotropy traces the large-scale structure (i.e., $). Various secondary 
effects have been studied in the liter ature: the Sunyaev-Zel'dovich effect [12|, cosmic 
reionization 124J, point sources 1251 ] . and ISW 121 ] . It has been shown that the last 
one, the ISW-lensing coupling, is the most dominant contamination of f^i? 1 

Using equation [46] in equation HU one finds S\ ens -s as 



(47) 
(48) 



14 



Slens-S =\J d 2 h{P(h)[d 2 E](h)Q(fl) 



- [d 2 P](fi)E(h)Q(h) - P(h)E(h)[d 2 Q]{h) 
+ (linear terms)} , 

which can be measured from the data. Here, the "linear terms" contain 9 terms with 
()mc, such as — P(h)([d 2 E](h)Q(h)) MC , etc. The map E{n) is given by equation |2T| 
and the other maps are defined by 

iPtn-x. 



P(n) = ^2wiCfei(h) =J2™i C f( C L a)imXim[,n) 



in i 



Q{h) 



J2uJiCf b ei (h) = Y.wiCr(C- l a) lm Y lm (n) 



(49) 
(50) 



Im 
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The maps with d 2 are given by d 2 P = — J2i + ^-) w iCf ei(n), etc. The map P(n) is a 
Wiener- filtered map of the primary temperature anisotropy from z = 1090. 



4-2. Lensing-ISW coupling 

A change in th e cur vature perturbation yields a secondary temperature anisotropy via 
the ISW effect [ToSj : 

' ^ dr ^(r, nr), (51) 



T Jo dr 

where r is the comoving distance and r* is the comoving distance out to z = 1090. Here, 
note again $ = — $Ncwton- The cross-power spectrum of and the ISW effect is then 
given by 

C*^ = trdr T -^P w [ l -A, (52) 
Jo r*r A \r J 

where P$$/(fc,r) is the cross-power spectrum of $ and $' = <9$/<9r, which can be 
calc ulat ed from the power spectrum of $, P$(/c,r), as i"W(A;,r) = |[9P$(fc, r)/9r] 



13 



126 



Here, P$(k, r) is not the primordial power spectrum, but it includes the linear 
transfer function, T(k), and the growth factor of $, g(r): 

P*(Kr) = j^[T(k)g(r)} 2 . (53) 

Using this, one finds P$$/(fc, r) = (g' / g)P$(k, r). Note that g>(r) is normalized such that 
g(r) = 1 during the matter-dominated era. 

With this result, it is easy to see why the lensing-ISW coupling yields the squeezed 
configuration: on very large scales, where T{k) — > 1, cf ,lsw oc l// 3 . On smaller scales, 
T{k) declines with k, and thus cf' lsw falls faster than 1/Z 3 . The lensing coupling 
includes 1(1 + l)^' 1 w , which falls faster than 1/Z, i.e., the largest power comes from 
the smallest I. 

A recent estimate by Hanson et al. [Tsj] showed that the lensing-ISW coupling, 
if not included in the parameter estimation, would bias /nl 8,1 by A/^L al = 9-3- The 
expected bias for WMAP is A/^ al = 2.7 [11]. One can remove this bias by including 
the lensing-ISW coupling (or any other lensing- secondary couplings) using the optimal 
estimator given by equation |48j 



5. Second-order effect 

5.1. General discussion 

So far, we have assumed that one can use equation [281 

/d^k 
^33$ P (fc)0rj(A;)y^(fc), 

to convert the primordial curvature perturbation to the temperature anisotropy. (Here, 
the subscript "p" stands for "primordial? by which we mean $ p = |£ without the 
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linear transfer function.) However, this equation is valid only for linear theory. As any 
non-linear effects can produce non-Gaussianity, one has to study the impacts of various 
non-linear effects on the observed non-Gaussianity. 

The origin of the linear radiation transfer function is the linearized Boltzmann 
equation: 

— — + ikfiA w + <x T n e aA (1) = S {1) (k, a, n), (54) 
or} 

where rj is the conformal time, \i = k ■ h, A^ = 4[ATW(fc, j u, v)/T] is the perturbation 
in the photon energy density, and is the linear source function, which depends 
on the metric perturbations as well as on the density, velocity, pressure, and stress 
perturbations of matter and radiation in the universe and the photon polarization. 
The second-order Boltzmann equation is then similarly written as 
,9A( 2 ) 

— — + ikfiA^ + a T n e aA^ = S {2) (k, n, 77), (55) 
or} 

where A (2) = 8[AT (2) (fc, n, r})/T] + 12 [AT (1) (k, fi, r/)/T] 2 , and S {2) is the second-order 
source function. Note that the azimuthal symmetry is lost at the second order, and thus 
the perturbations depend on the directions of k and h independently. In this case, the 



second-order ai m is given by [127 



X 



Y,Ft\k',k",k)Y* ml {k), (56) 



where F}^ 1 is the second-order radiation transfer function, whose form is determined 
by the second-order source function, S^ 2 \ in the Boltzmann equation. 

The shape of the second-order bispectrum, {a\l m2 a\l m3 a\l m3 } , is determined by the 
shape of the second-order radiation transfer function. If the second-order radiation 
transfer function vanishes in the squeezed limit, i.e., F}^ 1 ' \k' , k" , k) — > for k — > 0, 
then the CMB bispectrum would not peak at the squeezed configuration, and thus the 
resulting /^l 3,1 would be small. 



The second-order source function is quite complicated 128Ml38j . but it can be 
divided into two parta^: 

(1) The terms given by the products of the first-order perturbations, such as [$^)] 2 . 

(2) The terms given by the "intrinsically second-order terms," such as $^ 2 \ 

The intrinsically second-order terms are sourced by products of the first-order 
perturbations, and thus it is created by the late-time evolution of cosmological 
perturbations, whereas the terms in (1) are set by the initial conditions. 

% This decomposition is not gauge invariant, and thus which terms belong to (1) or (2) depends on the 
gauge that one chooses. Therefore, one must specify the gauge when making such a decomposi tion. Our 
discussion in this section is based o n th e gau ge choice made by Bartolo, Matarrese and Riotto 128l 129j 



and Pitrou, Uzan and Bernardeau 1361 . Il37j | , which reduces to the Newtonian gauge at the lin ear o rder. 
This seems a convenient gauge, as the products of the first-order terms only give I/nl^I < 1 127]. 
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The contamination of /nl^ due t o the terms in (1) is small, I/nl^I < 1 127 . 



Recently, Pitrou, Bernardeau and Uzan 137j reported a surprising result that the terms 
in (2) would give f^ 1 ~ 5 for the Planck data (7 max = 2000). 

Why surprising? As the intrinsically second-order terms arise as a consequence 
of the late-time evolution of the cosmological perturbations, they are generated by 
the causal mechanism, i.e., gravity and hydrodynamics. It is difficult for the causal 
mechanism to generate the bispectrum in the squeezed configuration, as it requires very 
long wavelength perturbations to be coupled to short wavelength ones. 

5.2. Newtonian calculation 

As an example, let us consider the well-known second-order solution for <3>( 2 ) in the 
sub-horizon limit, i.e., k 3> aH, which is equivalent to taking the non-relativistic 
(Newtonian) limit. Here, the second-order Barde en curvature perturbation is defined 
by $ = $(!) + i$( 2 ). The explicit solution is 



136 



k'k"r)^ 
k 



x F^ ) (k',k")^ 1 \k')^ 1 \k"), (57) 
where the linear perturbation, $W, on the right hand side is constant during the matter 



(s) 

dominated era, and the symmetrized function, F 2 , is defined as 

<'), 




2 



Fr(k u k 2 ) = ^+--^-^ \ ^ + ^ +=; ■ (58) 



as 



Note that F 2 is related to the function G given in equation (8.9) of 129 
G(k x ,k 2 , k) = (^) 2 F 2 w (fci, k 2 ) with k = fci + k 2 . 

The function F 2 s \ki,k 2 ) vanishes in the squeezed limit, ki = —k 2} and thus the 
CMB bispectrum generated from $( 2 ) in the Newtonian limit is not given by the local 
form. To see this, let us calculate 

($(&!, 77)$(fc 2 , vMk 3 , rj)) = (2Tr) 3 5 D (k 1 + k 2 + k 3 )F 2nd (k 1: k 2 , fc 3 , rj), (59) 



where 



rj 2 

F 2n d(ki, k 2 , k 3 ,rj) = — 



kik 2 s 



F 2 (s) (fc 1 ,fc 2 )P $ (fc 1 )P $ (fc 2 ) + (2 perm.) 



(60) 



and P$(fc) = AT 2 (k)/k A - n *. The shape dependence of F 2nd (fci, k 2 , fcsX^AO^Ai) 2 is 
shown in Figure H]for various values of k\ (because P 2n d is not scale invariant). The CMB 
data are sensitive to k\ < A; max ~ 0.2 h Mpc _1 (/ max /2000). We find that the bispectrum 
peaks at the equilateral configuration on large scales {k\ < 10~ 2 h Mpc -1 ), and it peaks 
along the elongated configurations on small scales (k% ~ 0.1 h Mpc -1 ). It peaks at 
the squeezed configuration on a very small scale (k\ ~ 1 h Mpc -1 ), but these scales 
are not accessible by the CMB due to the Silk damping. Note that the most squeezed 

+ Note that our $ is (—1) times $ used in equation 72 of [136]. 
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configuration shown in this Figure has k% = k-i — IOO/C3. The dominant shape changes 
with scales, as the linear transfer function, T(k), declines with k, with the small-scale 
limit given by T(k) oc \nk/k 2 . From these results, we expect the second-order effect in 
the Newtonian limit to yield only a small contamination of /^l ■ 

The dominant contribution to the second-order temperat ure a nisotropy in the sub- 
horizon limit is given by the second-order Sachs- Wolfe effect 



136 



AT® 



n 



I ■ ■ - 1 .) ■ ■ ' - \ • " ■ ■ • ' " 1 ■ (61) 

where -R* = 3pf,/(4p 7 ) is the baryon-photon ratio at the decoupling epoch. (Here, a 
factor of 1/2 comes from our way of defining the second-order temperature anisotropy, 
AT = AT^ + AT^ and the se cond -order curvature perturbation, $ = $W + |$( 2 ). 
This definition follows from Ref. 1271 ].) The corresponding second-order a\ m is 




d 3 k"5 D (k' + fe" - fe)$W(fc / )^ 1) (fe") 



{k'X)nti)T{k")j l {ku)8 lv 5 n 



(62) 

where the linear primordial perturbation, is related to as $W(fc) = 

&^\k)T(k). Comparing this with equation [56l we identify the term inside the square 
bracket as the second-order radiation transfer function, F/^ 1 '^] 

With this result, one can calculate the reduced bispectrum of the Newtonian second- 
order effect, bf^ 3 . The resulting f^ 3,1 is always less tha n un i ty re gardless of the angular 



scales (D. Nitta 2010, private communication; also see 140l . Il4l| ). The calculation was 
done for I < 2000. 



How can we reconcile this result with those found in 1371 ] ? The calculations given 
above (equations |57] and 1ST]) are valid only in the sub-horizon limits, and thus they are 
not suitable for calculating the contribution to the squeezed-limit bispectrum, which can 
correlate super- and sub-horizon fluctuations. Therefore, the difference between these 
results seems to imply: 

1. The dominant contamination of /^l comes from the general relativistic (post 
Newtonian) evolution of $ < - 2 - ) that is not captured by the above Newtonian 
calculation (equation 1571). 

2. The full second-order radiation transfer function beyond the sub- horizon 



approximation (equation |6TJ) gives the dominant contribution to f- 



local 
NL • 



The CMB bispectrum from the second-order ISW effect was considered in 139, I I 



Incidentally, in the notation of 
\R* (^^) 2 F 2 (s) (fc',fe")T(fc')T(fc")- 



[l27l| (see their equation 2.29), 



Sqo (^'i V*) 
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Figure 4. Shapes of the second-order bispectrum due to the second-order curvature 
perturbations in the Newtonian limit given in equation I60I Each panel shows the 
normalized amplitude of i*2nd.(&i) k2, k^)(k2 / ki) 2 {ky, / ki) 2 as a function of fca/fei and 
k^/ki for a given k\, with a condition that k% < k2 < k\ is satisfied. The amplitude 
is normalized such that it is unity at the point where -F2nd(^i, &2 ; k^){k2 / ki) 2 {k^ / k\) 2 
takes on the maximum value. (Top Left) k x = 10 -3 h Mpc -1 . (Top Right) h = 
IQ- 2 h Mpc -1 . (Bottom Left) ki = 10 -1 h Mpc -1 . (Bottom Right) k\ = 1 h Mpc -1 . 
The CMB data are sensitive to k\ < fc max ~ 0.2 h Mpc -1 (/ max /2000), where the 
second-order bispectrum peaks at the equilateral configuration on large scales, and 
peaks along the elongated configurations on a smaller scale {k\ ~ 0.1 h Mpc -1 ). On a 
very small scale (fci ~ 1 h Mpc -1 ), it peaks at the squeezed configuration. Note that 
the most squeezed configuration shown here has k\ = k% = 100fc3. 



Perhaps both contributions are important. This is yet to be confirmed; however, if this is 
true, one should be able to construct a simple template for the second-order bispectrum, 
and use it to remove the contamination by including its amplitude, in the fit. 

6. Four-point function: local-form trispectrum test of multi-field models 

Widely used notation for the "local-form trispectrum" is 

($(fe 1 )$(fc 2 )$(fe 3 )^(fc4)> 

= (2n) 3 5 D (k 1 + k 2 + k 3 + k 4 ) 

x {^r NL [P*(fci)P*(fc 2 ) {P*(fc 13 ) + P$(ku)} + (11 perm.)] 

+ 60nl [P$(ki)P${k 2 )P*(h) + (3 perm.)]}, (63) 

where k^ = |fcj + kj\. (In this section, we use $ for the primordial perturbation, i.e., 
$ = $ p = When the curvature perturbation is given by the simplest local form, 
$ = $ L + /nl 11 *! + 9NL$h one fi nc is the above trispectrum with t nl = (6f^l al /5) 2 



44l | . However, in general tnl is different from (6/^l /5) . 
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To see this, let us consider a broad class of multi-field models in which the primordial 
curvature perturbation, (, is given in terms of the field derivatives of the number of e- 
folds, iV = In a, and the perturbation in the J-th scalar field, 6<pi, as 



c = E 



ON 



z IJ 



o 2 n 



ij + • • • 



where (S^jSSj) 



for / ^ J 



(64) 



This expansion is known as the "5N formalism" 



22, M, ll42|, ll43|. In this case, and r NL are given by [22], M, ll44 



6 r local 



' NL 

E/(iV,/) 2 ] 2 ' 
EijkNijNjNjkNk 



1213 



^213 



where Nj = dN/d<p! and N u = d 2 N/d(f) I d(f) u 



[J2i(Nij _ 
Suyama and Yamaguchi 



49 



that the Cauchy-Schwarz inequality implies that the following inequality, 



TNL > 



is satisfied. To derive this result, use the Cauchy-Schwarz inequality: 




with 



aj 



[EANjf]^' 
[Ej(^,j) 2 ] 1/2 ' 



(65) 

(66) 
showed 

(67) 

(68) 

(69) 
(70) 



is satisfied for the simplest local-form model, ( 



The equality, r NL = (6/^° L cal /5) 2 
NxSfa + (N 11 /2)5<Pi 

Note, however, that one finds a different relation between and t nl when the 
Cauchy-Schwarz inequality becomes trivial, i.e., = 0. For example, when ( is given by 

18 2 N 



44 



c 



dN 



+ 



(71) 



0, one finds tnl ~ 10 3 (/NL al ) 4/ ' 3 [50|] . The Cauchy-Schwarz inequality 



uipi 2 

and (5(f) l&fo) 

becomes = because aj = for all /, and thus both the bispectrum and the 
trispectrum come from the second term in equation [7U the bispectrum is given by 
the 6-point function of 5<p2, and the trispectrum is given by the 8-point function of 5<p2- 
In this case, whether tnl > (6/NL al /5) 2 is satisfied depends on the value of f^ al . For 
this particular example, the current limit of /j^l < 74 implies that tnl > (6/NL al /5) 2 
is still satisfied parametrically. 



The inequalit y i s valid also for the "quasi-single field inflation" of 145 



equation (7.6) of [69|. (Note that rf^ and t^ 1 l in 
in this paper, respectively.) 



sec 



69[ correspond to tnl and g^L 
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Bciiid 


TTnrp arm inn 


£ local 
J NL 


i-cquil 

Jnl 


i-ortiiog 

Jnl 


h 

u src 


v+w 


Raw 


59 ±21 


33 ± 140 


-199 ± 104 


N/A 


v+w 


Clean 


42 ±21 


29 ± 140 


-198 ± 104 


N/A 


v+w 


Marg. 


32 ±21 


26 ± 140 


-202 ± 104 


-0.08 ±0.12 


V 


Marg. 


43 ±24 


64 ± 150 


-98 ±115 


0.32 ±0.23 


w 


Marg. 


39 ±24 


36 ± 154 


-257 ± 117 


-0.13 ±0.19 



Table 1. Estimates and the corresponding 68% intervals of the primordial non- 
Gaussianity parameters (f^ al , /nl"' > /nl ° g ) an d *he point source bispectrum 
amplitude, 6 src (in units of 10 -5 /iK 3 sr 2 ), from the WMAP 7- year temperature maps. 
This table is adopted from 



Therefore, if observations indicate tnl < (6/NL al /5) 2 , a broad class of multi-field 
models satisfying the above conditions would be ruled out. This property makes the 
trispectrum a powerful probe of the physics of multi-field models. Note that the simplest 
local-form limit, tnl = (6/NL al /5) 2 j has no special significance, as this is just one of many 
possibilities of multi-field models. (All single-field models predict a non-detectable level 
of the primordial f^jf 1 , and thus the observational test using the above relation between 
/nl 3,1 an d t nl has no relevance to single-field models.) 

The expected 95% uncertainties in tnl from the 7-year WMAP data (7 ma ,y ~ 500) 
and the 1-year Planck data (Z max ~ 1500) are 5000 and 560, respectively |73|. // the 
Planck finds f^i? 1 ~ 30, then it would be able to test if the measured tnl would satisfy 
equation [#7| This provides an excellent science case for the trispectrum that would be 
measured by Planck. 

The expected uncertainties in #nl have not been calculated yet, although we expect 
them to be much greater than those for r NL , as #nl is the co effic ient of the cubic-order 
term (i.e., #nl is much more difficult to constrain than tnl) [146 . 



The local-form trispectrum is not the only possibility. Various other inflation 
models predict distinctly different quadrilateral shape dependence. For some general 



analyses of shapes, see |35|, |42|, |58|, |69 



Finally, while we do not discuss the large-scale structure of the universe in this 
article, the most promisin g probe of the local- form trispectrum seems to be the 



bispectrum of galaxies. See [86|, Il47l . Il48j for details. 



7. Current Results 



7.1. Bispectrum 



(Most of this subsection is adopted from Section 6.2 of [111].) 



In 2002, the first limit on /^l was obtained from the COBE 4-year data [1491 ] by 



150], using the angular bispectrum. The limit was improved by an or der of magnitude 
when the WMAP first year data were used to constrain /nl &1 151] . Since then the 
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limits have improved steadily as WMAP collects more years of da ta and the bispectrum 



152-156 



method for estimating /n^ 3,1 ^ as improved 

Using the optimal estimators described in Section 3, we have constrained the 
primordial non-Gaussianity parameters as well as the point-source bispectrum using th e 
WMAP 7-year data. The 7-year data and results are described in Refs. 



We use the V- and W-band maps at the HEALPix resolution = 1024. As the 
optimal estimator weights the data optimally at all multipoles, we no longer need to 
choose the maximum multipole used in the analysis, i.e., we use all the data. We use 
both the raw maps (before cleaning foreground) and foreground-reduced (clean) maps to 
quantify the foreground contamination of /nl parameters. For all cases, we find the best 
limits on /nl parameters by combining the V- and W-band m aps, and marginalizing 



over the synchrotron, free-free, a nd d ust foreground templates 1581 ]. As for the mask 



we always use the KQ75y7 mask [158 



In Table 17.11 we summarize our results: 

1. Local form results. The 7- year best estimate of /nl i s 

/nT 1 = 32 ±21 (68% CL). 

The 95% limit is —10 < /nl < 74. When the raw maps are used, we find 
/nl 11 = 59 ±21 (68% CL). When the clean maps are used, but foreground templates 
are not marginalized over, we find /nl 3,1 = 42 ± 21 (68% CL). These results (in 
particular the clean-map versus the foreground marginalized) indicate that the 
foreground emission makes a difference at the level of A/^^ 1 ~ 100 We find 
that the V±W result is lower than the V-band or W-band results. This is possible, 
as the V±W result contains contributions from the cross-correlations of V and W 
such as (WW) and (VWW). 

2. Equilateral form results. The 7-year best estimate of /nl" 1 is 

/ N q L uil = 26 ± 140 (68% CL). 

The 95% limit is -214 < / N l uil < 266. For f^ 31 , the foreground marginalization 
does not shift the central values very much, A/ N l Ui1 = —3. This makes sense, as 
the equilateral bispectrum does not couple small-scale modes to very large-scale 
modes I < 10, which are sensitive to the foreground emission. On the other hand, 
the local form bispectrum is dominated by the squeezed triangles, which do couple 
large and small scales modes. 

3. Orthogonal form results. The 7-year best estimate of /nl^ 8 is 

^orthog = _ 2Q2 ± (6g% 

The 95% limit is —410 < og < 6. The foreground marginalization has little 
effect, A/°f os = -4. 

ff The effect of the foreground marginalization depends on an estimator. Using the needlet bispectrum, 
Cabella et al. [l6l| found /^°£ al = 35 ± 42 and 38 ± 47 (68% CL) with and without the foreground 



marginalization, respectively. 
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As for the point-source bispectrum, we do not detect 6 src in V, W, or V+ W. In [85[ , 
we estimated that the residual sources could bias by a sma U positive amount, and 
applied corrections using Monte Carlo simulations. In this paper, we do not attempt to 
make such corrections, but we note that sources could give A/^° L cal ~ 2 (note that the 
simulations used by [85] likely overestimated the effect of sources by a factor of two). 

extrapolating the previous results 



85 



85 



but we have 



As the estimator has changed from that used by 
is not trivial. Source corrections to f^ 11 and j^J" 6 could be larger 
not estimated the magnitude of the effect for the 7-year data. 

As we described in Section 4, among various sources of secondary non-Gaussianities 
which might contaminate measurements of primordial non-Gaussianity (in particular 
/nl^) > a coupling between the ISW effect and the weak gravitational lensing is the most 
dominant source of c onfu sion for /nl^- Calabrese et al. jl62| used the skewness power 
spectrum method of 1631 ] to search for this term in the WMAP 5-year data and found 
a null result. 



7. 2. Trispectrum 

The optimal estimators for the trispectrum have not been implemented, largely because 
they are computationally demanding. While the first mea surem ents of the angular 
bispectrum were made from the the COBE 4-year data 149 ] by Il64 | in 2001, limits 
on the physical parameters have not been obtained from the direct trispectr um a nalysis. 

Recently, Smidt et al. 7_5] used the sub-optimal estimator developed in [1651 ] (which 
becomes optimal in the limit that the instrumental noise is isotropic) and found the 
95% CL limits of -7.4 x 10 5 < g NL < 8.2 x 10 5 and -0.6 x 10 4 < r NL < 3.3 x 10 4 . The 
current limit is consistent with the Suyama-Yamaguchi inequality, tnl > (6/NL al /5) 2 . 



7.3. Other statistical methods 

While the optimal estimators for the /nl parameters (with the minimum variance) 
must be constructed from the PDF as in Section 3, there are various other ways of 
constraining non-Gaussianity. While these other methods are usually sub-optimal, they 
serve as useful diagnosis tools of the results obtained from the direct bispectrum and 
trispectrum methods. In some cases, they are easier to implement than the optimal 
estimators. 



A major 



jrogr ess in the topological Gaussianity test using the Minkowski 



junctionals 166l4l70| since 2004 is the derivation an d impleme ntation of the analytical 



formula for the Minkowski functionals of the CMB 171. 172 . This method has been 



applied to the WMAP data ll73j, ll74j as well as to the BOOMERanG data [175j]. The 



local 



Planck data are expected to reach the 68% limit of A f NI 



20 1711, which is worse than 



82 



An advantage of the Minkowski 



the limit from the optimal method, Af^l = 5 
functionals is that the measurements of the Minkowski functionals do not depend on 
the models, and thus the computational cost is the same for all models. This allows 
one to obtain limits on various models, for which the optimal estimators are difficult to 
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implement. For example, a limit on the primordial non-Gaussianity in t he iso curvature 



174 



perturbation is currently available only from the Minkowski functionals 

Instead of expanding the temperature anisotropy into spherical harmonics, one 
may choose to expand it using a different basis. One popular basis used in the C MB 
community is the so-called Spherical Mexican Hat Wavelet (SMHW). See |176l . Il77| for 
reviews on this method. A major progress in this method is the realization that the 
3-point function of the wavelet coefficients made of large and small smoothing scales is 
nearly an optimal estimato r for the local- form bispectrum: when only the adjacent scales 
are included, Curto et al. H found -8 < /^ al < 111 (95% CL) from the WMAP 5- 
year data. When all the scales (including large-small scale combinations) are incl uded in 
the analysis, the limit improved significantly to — 18 < f^ 8,1 < 80 (95% CL) 1791 ]. which 
is similar to the optimal limit from the 5-year data, —4 < /nl^ < 80 (95% CL) [10 . 
An advantage of the SMHW is that it retains information on the spatial distribution of 
the signa l. This property can be used to measure f^l al as a function of positions on the 
sky jl79| . In addition, the analytical formula for the SMHW as a function of f^ al has 



180 



181 



been derived and implemented (A. Curto 2009, private communication). See 
for earlier limits on f^i? 1 from the SMHW. 

Another for m of spherical wavelets that has been us ed to constra in f^L* 1 is the 
spherical needlets 1821 ] . The limits on /^ al are reported in 16ll . ll83l . ll84| . This method 
also allow s one to look for a spatial variation in f^f^ 1 , and the results are re ported in 



185 



186 



For the other types of wavelets considered in the literature, see 



187. 



188 



and references therein. 

Many other statistical methods have been propo sed and u sed for constraining /nl 
in the lite ratu re. An incomplete list of referen ces i s: 1891 . Il90| on the real-space 3- point 
function; 1911 ] on the integrated bisp ectrum; 1921 ] on the 2-1 cumulant correlator; 181 



iunction; on tne integrated bisp ectrum; [l^zl] on tne z-1 cumulant correlator; [ioi 
on the local curvature; and 0, H on the iV-point PDF. Also see references therein 
While we have not listed the statistical methods that have not been used to constrain 
the primordial non-Gaussianity parameters yet, there are many other methods proposed 
in a general context in the literature. 



8. Conclusion 



Since the last review articles on signatures of primordial non-Gaussianity in the CMB 
were written in 2001 [l| and 2004 [2], a lot of progress has been made in this field. The 
current standard lore may be summarized as follows: 

1. Shape and physics. Different aspects of the physics of the primordial universe 
appear in different shapes of three- and four-point functions. 

2. Importance of local shape. Of these shapes, the local shapes have special 
significance: a significant detection of the local-form bispectrum (with f^ &l ^ -0 
would rule out all single-field inflation models, and the local-form trispectrum can 
be used to rule out a broad class of multi-field models (if not all multi-field models) 
by testing r NL > (6^° L cal /5) 2 . 
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3. Optimal estimators. The optimal estimators of the bispectrum and trispectrum 
can be derived systematically from the expansion of the PDF. The optimal 
bispectrum estimator has been implemented. 

4. Secondary. The most serious contamination of /^l * s due t° the lensing-ISW 
coupling, which can be removed by using the template given in Section 4. 

5. Foreground. The Galactic foreground contamination is minimal for /nl" 1 an d 
/NL h0g 5 Di rt it can be as large as /nl^ ~ 10 for the local- form bispectrum. This 
must be carefully studied and eliminated in the Planck data analysis. The random 
(Poisson) point-source contamination can be removed by using the template given 
in Section 3.2. 

Some outstanding issues for the "CMB and primordial non-Gaussianity" include: 

1. Second order. (In Newtonian gauge) the products of the first-order terms and 
the intrinsically second-order terms in the sub-horizon limit do not contaminate 
the local-form bispectrum very much (A fj,° L cal < T). However, would the post- 



Newtonian effect give A/^L al ~ 5, as found by [1371 ]? If so, we need to construct a 
template for this effect. 

2. More foreground. How can we model the non-Poisson (clustered) point source 
bispectrum? How about the foreground (and secondary) contamination of the 
primordial trispectrum? 

3. Trispectrum estimators. How can we implement the optimal trispectrum 
estimators for both local and non-local shapes? 

These issues would become important when the Planck data are analyzed in search of 
primordial non-Gaussianity. The Planck is expected to reduce the uncertainty in f^i^ 1 
by a factor of four compared to the current limit, f^ al = 32 ± 21 (68% CL). If the 
Planck detected /^L al ~ 30, then the trispectrum would provide an important test of 
multi-field models. In particular, if f^l al ^ 1 aric ^ t nl < (6/NL al /5) 2 are found, then 
all single-field models and many (if not all) multi-field models would be ruled out, and 
thus the standard paradigm of inflation as the origin of fluctuations would face a serious 
challenge. 

However, do not despair even if the Planck did not detect the primordial bispectrum 
or trispectrum - while the CMB may end its leading role as a probe of primordial 
non-Gaussianity (unless the next-generation, comprehensive CMB satellite which can 
measure b oth the temperature and polarization to the cosmic-variance-limited precision 



is funded [l94j ). the large-scale structure of the universe would eventually take over and 
substantially reduce the uncertainties in the local-form parameters such as f^l 31 , t nl , 
and #nl (see Desjacques and Seljak's article in this volume). 
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